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Introduction. Transverse vibrations of a bimorph consisting of two piezomagnetoelectric layers and located in the 
alternating magnetic field are investigated. Piezomagnetoelectric layers are multilayer composites with alternating 
piezoelectric and piezomagnetic layers. The mechanical and physical properties of such a composite are given by 
known effective constants. 

Materials and Methods. The applied theory of multilayer plate vibrations takes into account the nonlinear distribution 
of electric and magnetic potential in piezoactive layers in the longitudinal and transverse directions. On the basis of this 
theory, the stress-strain state, the dependences of deflection, electric and magnetic potentials on the volume ratio of the 
composition of the hinged bimorph, are investigated. The electric potential 1s assumed to be zero at all electrodes, while 
the magnetic potential is zero at the inner boundary and unknown at the outer boundaries. Therefore, the distribution of 
electric and magnetic potentials in the middle of the layer are unknown functions. In the case of the magnetic potential, 
the distribution at the outer boundary is also unknown. In the problem, the Kirchhoff hypotheses for mechanical 
characteristics were accepted. The use of the variational principle and the quadratic dependence of the electric and 
magnetic potentials on the thickness of piezoactive layers made it possible to obtain a system of differential equations 
and boundary conditions. 

Results. When the volume ratio of the composition of piezoactive bimorph materials changes, the electric potential in 
the middle of the layer changes nonlinearly. The magnetic potential in the middle of the layer and at the outer boundary 
increases almost linearly with an increase in the volume percentage of Ba7iO3. The dependence of the deflection in the 
middle of the layer is determined. 

Discussion and Conclusions. An applied theory for calculating transverse vibrations of a bimorph with two 
piezomagnetoelectric layers is constructed. The dependence of the characteristics of the stress-strain state, electric and 


magnetic fields on the volume fractions of piezomagnetic and piezoelectric materials, is investigated. 
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Introduction. In the production of sensor and measuring systems, small household appliances, cell phones, 





and wireless sensor systems, powerful energy sources are not needed to monitor and diagnose the technical condition of 
objects. In this case, the prerequisites are the mobility and volatility of the above devices. 

Piezoelectric materials directly convert electrical energy into mechanical energy and vice versa. This property 
allows them to be widely used in science and technology. These materials are used in ultrasonic emitters of elastic and 
acoustic waves, receivers of such waves, devices for suppressing vibrations of machine elements and structures, etc. 
Recently, another field of application of piezoelectrics has been rapidly developing — energy collection and storage 
devices. In this case, piezoelectric materials are part of piezoelectric energy generators (PEG). PEG are placed on 
elements of machines or structures that vibrate intensively, which are in the zone of elastic wave propagation or are 
exposed to variable pressure. The main types of these devices have a bimorph or stack multilayer structure and 
experience bending or longitudinal deformations, respectively. Low-power sources of electric current are created on the 
basis of PEG. They include autonomous power sources (e.g., for damage monitoring devices in hard-to-reach places of 
pipeline structures, etc.). An overview of such devices is available in [1—2]. One of the ways to design effective PEG is 
the use of piezoactive composites of various types of connectivity and heterogeneous materials based on piezoceramics, 
including porous one. 

PEG, in whose design there are additional electromagnetic elements or permanent magnets, can fix or use the 
energy of an alternating magnetic field. One of the ways to solve this problem is the use of piezomagnetic materials in 
combination with piezoelectric ones. In this case, the alternating magnetic field causes the deformation of the 
piezomagnetic and the coupled piezoelectric, as a result, the latter generates electrical energy. There is a class of 
materials with ferromagnetic properties. Piezomagnetism is a phenomenon observed in some antiferromagnetic and 
ferromagnetic crystals. It is characterized by a linear relationship between the magnetic polarization of the system and 
mechanical deformation. In a piezomagnetic material, a spontaneous magnetic moment can be induced through 
applying mechanical stress, or deformation by applying a magnetic field. In studies of piezomagnetic materials [3-5], 
CoFe,O41s very often considered. In [6-8], a composite based on CoFe,O, and BaTiO; with piezoelectric and 
piezomagnetic properties is studied simultaneously. 

Solutions to the problems of electroelasticity and magnetoelasticity are given in [9-11]. In [12], applied 
theories of vibrations of multilayer piezoelectric plates were developed considering the specifics of the distribution of 
electric potential over the structure thickness. 

The problems on steady-state oscillations of an electro-magneto-elastic layer and a half-space under the action 
of harmonic loads are presented in [13, 14]. The prestressing is taken into account, as well as various electrical and 
magnetic conditions at the boundaries. The effect of these factors on the dispersion properties is investigated. 

Earlier [15, 16], an applied theory was developed that considers the inhomogeneous distribution of the electric 
potential in the longitudinal direction, and the quadratic dependence on thickness. In the same papers, the stress-strain 
and electrical state of a hinged and cantilevered bimorph is investigated. In both cases, the applied theory showed good 


convergence with the finite element modeling results. The authors also developed an applied theory of bimorph 
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vibrations [17] consisting of an electroelastic and magnetoelastic layer. This approach is in good agreement with the 
finite element analysis results. 

In this paper, the vibrations of the device are considered in the framework of a plane deformation. Based on the 
variational principle, an applied theory of bending vibrations of a two-layer piezoelectric bimorph is constructed. For 
steady-state vibrations, boundary conditions and a system of differential equations are obtained for four unknown 
functions (deflection, electric potential in the middle of the layer, magnetic potential in the middle of the layer and at 
the outer boundary), depending on the length of the bimorph. The influence of different percentage volume ratios of the 
bimorph composition on deflection, electric and magnetic potentials in certain positions, is investigated. The study 
results provide selecting the makeup of a composite piezomagnetoelectric material to achieve the most efficient 
operation of the device. 

Materials and Methods. A plate consisting of two identical piezoelectric layers is considered. It performs 
steady-state transverse vibrations within a plane deformation. Each layer is a 2—2 connectivity composite consisting of 


alternating piezoelectric and piezomagnetic layers (Fig. 1). 








Piezoelectric layer =a Piezomagnetic layer 


Fig. 1. 2-2 composite structure 
Effective properties of such a composite were found in [8]. Large surfaces of the layers are electrodated, and 
the layers themselves are polarized in thickness. The bimorph is hinged at the edges, all surfaces are free from 
mechanical stresses. The upper and lower boundaries of the plate are affected by magnetic flux B,, while at the 
boundary between the layers, the magnetic potential 1s considered to be zero. The electrical potential is zero at all 
electrodes. The side surfaces are considered to be insulated from magnetic and electric fields. 
The equations for describing the vibrations of a composite with effective properties, connectivity of 
mechanical, electric and magnetic fields, have the form [18]: 
V-o+pf =pu, V-D=0,, V-B=0, 
G6=cit-e -B=ir -H, 
D=e:6+k-E+a-A, 
B=h:e+a' -E+u-H, 


e=—(Vu+(Vu)'), Hao, Fs. (1) 
Here, O and € — mechanical stress and strain tensors, D and E — vectors of electric induction and electric field 
strength, B and H — vectors of magnetic induction and magnetic field strength, p — material density, c — elastic 
moduli tensor, e — piezoelectric moduli tensor, / — piezomagnetic moduli tensor, « — dielectric permittivity 
tensor, o& — magnetoelectric moduli tensor, 4 — magnetic permeability tensor, / — mass force density vector, 6, 
— electric charge volume density, u — displacement vector, @ and € — electrical and magnetic potentials. 


The boundary conditions are determined for the mechanical, electric and magnetic fields, respectively. 


For the first case, we note the absence of mechanical stresses at the bimorph boundary: 


o,,-7,| 


The biomorph is hinged at the ends (Fig. 2): 


UU, f7 13: 
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Fig. 2. Geometry and boundary conditions of bimorph with composite piezomagnetoelectric layers 


Next, we formulate the electrical boundary conditions. Electrical potential on the internal and external 
electrode, respectively: 


| =f = Vi. Ol, i _ V, 7 


We indicate the magnetic boundary conditions. Magnetic potential at the inner boundary: 
é|  =M,. 


x3 =0 


Magnetic flux 8B, affects the upper and lower boundaries of the plate: 
=B 


I. =tH 0° 
We use the variational equation for steady-state vibrations [10]. It generalizes Hamilton's principle in the 
electroelasticity theory taking into account magnetic components. For the case of plane deformation in the absence of 


surface loads and in the presence of magnetic flux: 


| | SHdS — pw’ | | u,du.dS + | | (p8u, +,59 + B,dE)dS =0, (2) 
S S Ss 


where 5H = Ogee D,6E;- B,5H;. 
To construct an applied theory of vibrations, we will accept Kirchhoff's hypotheses. In accordance with them, 


the distribution of displacements along the thickness has the form: 


Uy (X1,X3) =-X3Wi1, Ug (X1,X3) = W(X) (3) 

In particular, the single normal hypothesis is accepted for a mechanical field. Next, we consider a problem in 

which the value of the electric potential on the electrodes can be zero, so its distribution is not described by a linear 
function. Taking into account the possible inhomogeneity in the length of the element associated with the influence of 


boundary conditions at the ends of the bimorph, its thickness distribution is assumed to be quadratic: 
Ma 4X; Ne 2X 
x,,X,;) =V,(x,) =| —-1]4+V, (x +V. <=} +41]. 4 
(x5) = so) [ 7 {- ls (x) =| 2 +1) (4) 


Here, X; = x3-h/2. Functions Vo, V, and V, are responsible for the value of the electric potential at the inner electrode, 
in the middle of the layer, and at the outer electrode, respectively. To satisfy the conditions of the problem, we take 
these functions in the following form (see Fig. 2): 
Vi(~%) =Y =const, V(x) =P), V(x) =V, = const. 
Here, function ®(x,) is unknown. 
We represent a quadratic distribution of the magnetic potential over the thickness of each layer. The 


distribution along the length is heterogeneous, at the inner boundary of the layers, its value is assumed to be zero: 


(x,,%,) = Mone [ ae eam cxp[- + Jeaseay3(22-1). (5) 
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Here, X3 = X3-—h/2. Functions My, M, and M, are responsible for the value of the magnetic potential at the inner 
boundary, in the middle of the layer, and at the outer boundary, respectively, and are taken as follows (Fig. 2): 
M,(~%) =M, =const, M,(x%,)==,(%), M,(~%) =2,(%). 
Here, functions £,(x,) and &3(x,) are unknown. 
We substitute relations (3)-(5) into equation (2) and integrate it by the bimorph thickness, and then we 
equate the coefficients with independent variations 5w, 6®, d=, and d#, to zero. Thus, we obtain a system of four 
differential equations (6) from four unknown functions depending on x, (then, we omit the subscript), and five 


boundary conditions (7). 
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* ‘ ‘ ‘ = o) ed ~ 
Here, the following designations were introduced: ¢, =G,—-¢3/C33, €31 = €31-C13€33/C33, h31 = 


ha,- C13h33/C33, O33 —— X%33- C33h33/C33, E33 = Con Oss" (Cae ° They occurred after satisfying condition 033 — 0 and 
exclusion of €33. 
Research Results. The results of the bimorph calculation according to the proposed theory are compared to the 


finite element calculation in the low-frequency region for the volume ratio of the piezoelectric and piezomagnetic 
components 80 % BaTiO; and 20 % CoFe 0,4. The comparison has shown that the error in finding the characteristics of 
the mechanical and magnetic fields is less than 1%. When determining the electric field in the middle part of the plate, 
the difference was about 5%. Describing the situation in the vicinity of the support points, it should be noted that the 
size of the neighborhood along the longitudinal coordinate 1s approximately equal to the bimorph thickness. Here, when 


determining the electric field, a difference of 20% is recorded. 
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The first step in studying the vibrations of a two-layer piezomagnetoelectric bimorph with a change in the 
volume ratio of BaTiO; and CoFe,O, in the composite is to determine its effective properties. Tables 1 and 2 present 
these properties found from the results of study [8]. 


Table 1 
Material constants (elastic modules, dielectric and magnetic permittivity) 
for different BaTiO; volume fraction 


Volume Elastic modules Dielectric permittivity | Magnetic permittivity 


fiction [GPa 10m TON STC 
coat ee eee ee 


0 | 2860 | 1730 | 1700 | 2695 | 45.30 | 0080 | 0.093 | 5.900 | 1570 
10} 2709 | 160.4 | 1549 [260.0 [45.07 | 1469 [0073 | 5318 | 0.632 


6 | 2067 | tos.t | 102.1 | 2156 | 4392 | 7490 | 0.220 | 2.390 | 0.155 
80 | 1860 | 923 | 859 | 2013 | 4346 | 9.488 | 048i 1200 | 0.120 
99 | 1766 | 854 | 789 | 1989 | 4323 | 10353 | 0.857 | 0635 | 0.110 


Table 2 





Material constants (piezoelectric, plezomagnetic and magnetoelectric modules) 
for different BaTiO; volume fraction 


Volume 
A A, SSS 
BaTiO; (%) eS SS eee ee 
of 
ea [ores ee es || 


x0 2a 


- | 0038 | 0220 [696 | 297 [228 | 20 | 34s | 2a 

80 | 0060 [057 [928 [137 91110 
[90 | 0263 | 1e7_| roa [46 [a6 ‘| 55 ([ 2 
[100 [4.400 | is600 [eo fo fo jo [0 | 0 





The bimorph vibrations were excited by a magnetic flux applied to the upper and lower faces (Fig. 2), which 
varied according to the harmonic law with an amplitude of B, =5x10~° Wb and a frequency of 10 kHz. 

Figure 3 shows the deflection in the middle of the layer depending on volume fraction of BaTiO. It can be 
seen from the graph that the deflection in the position having coordinates x, =L/2, x, =H/2, 1s zero if the bimorph 
consists only of a piezoelectric BaTiO;. The deflection of the bimorph reaches the greatest value if it contains only 


piezoelectric magnet CoFe,O,4. The deflection almost linearly depends on the volume ratio of the components of 


piezoactive materials. 
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Fig. 3. Deflection w(x, ) in the middle of the layer for different BaTiO; volume fraction 


Based on the data in Figure 4, it can be concluded that the electric potential in the middle of the layer varies 
nonlinearly with a change in the volume ratio of the composition of bimorph piezoactive materials. If the bimorph 
consists only of BaTiO; or CoFe,O,, then the electric potential at the point (L/2,H/2) is zero and reaches the highest 


value at 35 % BaTiO; in the bimorph. 
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Fig. 4. Electric potential ‘Y(x,) for different BaTiO; volume fraction 


The analysis of Figures 5 and 6 allows us to conclude that the magnetic potential in the middle of the layer =, (L/ 2) 


and at the outer boundary =,(L/2) increases almost linearly with an increase in the volume of BaTiO; in the bimorph. 
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Fig. 5. Magnetic potential &,(L/2) for different BaTiO; volume fraction 
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Fig. 6. Magnetic potential &,(L/2) for different BaTiO; volume fraction 


Discussion and Conclusions. An applied theory is proposed for calculating transverse vibrations of a bimorph 
made of two layers of a composite based on CoFe,O, and BaTiO; with both piezoelectric and piezomagnetic properties, 
in an alternating magnetic field. Such a design can serve as a model of a piezoelectric generator of a device for 
collecting and storing energy under the action of an external magnetic field. In the low-frequency region (below the 
natural frequency of the first bending mode), calculations of the stress-strain state of the bimorph, the distribution of 
electric and magnetic fields, are carried out. The dependence of the deflection, electric and magnetic potentials on the 
volume ratio of the bimorph composition is investigated. In further work, it 1s assumed to determine the output potential 
and power of an electric current excited by an alternating magnetic field. The purpose of these surveys will be to collect 
electrical energy. 
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